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In this study, recently introduced generalized distribution functions are 
summarized and by using one of these distribution functions, namely gener- 
alized Planck distribution, an alternative approach to the generalized Planck 
law for the blackbody radiation has been tackled. The expression obtained 
is compared with the expression given by C. Tsallis et al. [Phys. Rev. E52, 
(1995) 1447], and it is found that this approximate scheme provides bounds 
to the exact result, depending on the values of q index. 

I. INTRODUCTION 

Although the necessity of nonextensive statistical thermodynamics has been very clear 
for a long time in many physical systems such as three-dimensional self-gravitating astro- 
physical objects [1], black holes and superstrings [2], Levy- type random walks [3], vortex 
problem [4] and dark matter [5] where long-range microscopic interactions are present; an in- 
creasing tendency towards nonextensive formalisms keeps growing nowadays, along two lines: 
Quantum-Group-like approaches and Generalized Statistical Thermodynamics (GST). 

GST is recently introduced by C. Tsallis [6,7] and then the formalism not only has been 
applied to numerous concepts of statistical thermodynamics [8-28] , but also prosperous for 
the physical systems [29-33] where extensive Boltzmann-Gibbs statistics fails. The detailed 
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reviews of the general aspects of the formalism can be found in [34], and some investigations 
of the subject from the mathematical point of view are now available in [35]. 

The formalism has been based upon two axioms; namely (i) the entropy of the system 
is given by [6] 

S q = -k (!) 

where k is a conventional positive constant, pi is the probability of the system to be in a 
microstate, W is the total number of configurations of the system and q is a new parameter 
which is especially called the Tsallis g-index; (ii) the g-expectation value of an observable O 
is given by 

w 

(0) s = Eft^ • (2) 
1=1 

On the other hand, GST contains Boltzmann-Gibbs Statistics as a special case when q — 1. 
Clearly S q is extensive if and only if q = 1 and (1 — q) can be interpreted as the measure of 
the lack of extensivity of the system [24]. 

II. GENERALIZED DISTRIBUTION FUNCTIONS 

Amongst the systems handled in the frame of GST, the classical and quantum gases 
have been investigated for the first time by two of the authors of the present study and 
the results have been presented in two seperate manuscripts [13,14]. In the first paper, the 
fractal inspired entropies of the classical and quantum gases have been found by making 
use of the statistical weights and the generalized distribution functions have been obtained 
by introducing the constraints related to the statistical properties of the particles into the 
entropy [13]. In the second paper, the Tsallis entropy have been used for the same purpose 
and the same conclusions for the generalized distribution functions have been attained but 
in this case the constraints concerned with the statistical properties of the particles have 
been introduced in the calculation of the partition functions. However, the factorization 



2 



of the partition function which was used in this calculations has been the main disputable 
point of the approximation. We claim that this approximate procedure provides a bound 
to the exact results, depending on the values of q. In order to prove this, let us take a very 
simple physical system which could be in bounded or free states with the energy levels A 
and B. For this system let us take the inequality 



[1 + (1 - q)(A + B)\—« ^ [1 + (1 - q)A}— [1 + (1 - q)B]—« 



(3) 



It is easy to rewrite this in exponential form 

f 1 



cxp 



1-q 



\n[l + (l-q)(A + B)]\? 



expj^ln [1 + (1 - q)A] + In [1 + (1 - 



(4) 



After some algebra it is straightforward to find 

f 1 



cxp 



ln[l + (l-g)A+(l-g)£H^ 



cxp 



1-q 



In 



1-q 

l + (l-q)A+(l- q)B + (1 - qfAB 



(5) 



If the left-hand side and right-hand side of this inequality are called "exact" and "ap- 
proximation" respectively, then it follows that exact>approximation for q > 1 whereas 
exact < approximation for q < 1. It must be emphasized that this situation is valid if and 
only if the product A.B is positive. In photon case, which is the subject of this study, A 
and B, namely energy levels of the system are defined to be hv/kT that is always positive. 

In addition to this, an inequality can be derived for the partition function within the 
frame of (3). The grand partition function reads : 



£ 

{n} 



l-/3(l-g)£(e r -Ai) 



r=l 
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where {r,} are the elements of the grand canonical ensemble. For the photon case (/i — 0), 
the eq.(6) becomes 



z q = Y.[^-^-m { n } 

{n} 



(7) 



where Tt{ ri } = J2T=i e r is the Hamiltonian of the system. Employing the inequality (3), 
another inequality for the partition function can be written down as 

E [1-/5(1 -9) (ei + '-' + Ol^ 

in} 

Efl- ^(1 - 9) e i]^ •-- [1-/3(1- g)e m ]^ • (8) 
{n} 

Since the free energy is a monotonic function of Z q , namely, 

1 z x - q - 1 

^ = "^^7- • <») 

it is straightforward to write an inequality for the free energy. 

In the frame of the above simple mathematical analysis, it is evident that our approx- 
imation scheme provides a lower or upper bound to the exact results, depending on the q 
values. 

Within this approximation procedure the generalized distribution functions are given by 

« B = [i-(i-?Me r -/^F (io) 

K>r = — i — (ii) 

[l-(l-g)/3(e r - /U )]- + l 



"r>f * = : ; — rrn (12) 
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[l-(l- g )/3(e r - ^5=1-1 
where MB, FD and BE stand for Maxwell-Boltzmann, Fermi-Dirac and Bose-Einstein, 

respectively. In a very recent effort [36], quantum statistics has been studied by means of a 

kinetic approach and the distribution functions corresponding to the Tsallis probability are 

found to be the same with the eqs.(10)-(12). It is easy to show that the standard distribution 

functions are recovered in the q — > 1 limit. On the other hand, when \x is set equal to zero 

in eq.(12), one can obtain the generalized Planck distribution (PD) 

K>f D = 1 (13) 

which is going to be used here for making an approximation to the derivation of the general- 
ized Planck radiation law by following the well-known procedure available in many textbooks 
[37] of Statistical Physics. 



III. THE GENERALIZATION OF THE PLANCK RADIATION LAW 



Very recently, the generalization of the Planck radiation law is obtained by Tsallis et 
al. [24] in order to see whether present cosmic background radiation is (slightly) different 
from Planck radiation law due to long-range gravitational influence. In the course of their 
investigations, making use of the partition function given by Z q Z x jl — |(1 — q)f3 2 (7i 2 ) 1 | 
(where Z\ and {H 2 ) 1 stand for the values of these quantities in standard Boltzmann-Gibbs 
statistics) in the f3{l — q) — > limit, Tsallis et al. have achieved the generalization of the 
Planck law given in the following expression for q = 1 [24] 



D q (v)h 2 c 3 x 3 



7(1--*) 



1 + (1 - q)x 



l + e~ x x 1 + 3e" 



2 (1 -e~ x y 



(14) 



8n(kT) 3 e 

where D q {y) is the photon energy density per unit volume, v is the photon frequency and x = 
hv/kT. Then they applied this expression to the cosmic microwave background radiation to 
test for deviations from Planck radiation law and found a 95% confidence limit of \q — 1| < 
3.6 x 10~ 5 from the data obtained via Cosmic Background Explorer Satellite by Mather et al. 
[38]. In the present work we're trying to generalize the Planck law by using an approximate 
scheme which seems simpler and more general (not necessarily q = 1). 

The distribution of photons among the various quantum states with definite values of the 
momentum hk/2n and energies hu can be given by generalized Planck distribution (eq.(13)) 

K>r=- — . (15) 



9-1 



- 1 



[1-(1-9)&J 

On the other hand, the number of quantum states of photons with frequencies between v and 
v + dv is 8irVv 2 dis/c 3 (V being the volume of the photon gas). It is clear that multiplying 
this quantity by eq.(15), the number of photons in this frequency interval can be determined: 

8?rV v 2 dv 



dN = 



l-(l-g) 



hu 
kT 



(16) 
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Therefore the photon energy in this interval is given by 

8nhV v 3 dv 



dE = 



i-(i-g) 



hu 
kT 
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(17) 



- 1 



and finally the photon energy density per unit volume is 



„ . . 8nh v 3 , . 

D q {y) = - rT , (18) 



hu 
' kT 
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which generalizes the Planck's law. The comparison of this expression with that given by 
Tsallis et al. (eq.(5) of their paper) is illustrated in Figure. It is observed in the figure that, 
at low frequencies the two D q {y) plots completely fit into one another. On the other hand, 
towards the frequency values where the maxima of the curves occur, the plots diverge from 
one another by a certain amount. This is an expected result since the implication of our 
approximation scheme shows itself as an upper (lower) bound to the exact result when q > 1 
(q < 1). It is worthwhile to imply here that the bound is exactly same as the one that has 
appeared in ref. [39] where an application of the generalized distribution functions has been 
discussed (see Fig.l and Fig. 3 of this reference). 

Moreover, in the x = hu/kT « 1 case, we verify that eq.(18) becomes 

8n{kT) 3 x 3 
q[U) " h 2 c 3 x + !2^ x 2 + ( 2 -^) x s +~7. ' (19) 

which generalizes the Rayleigh- Jeans law. As it is expected, in the q — > 1 limit this expression 
transforms to D\(y) oc x 2 oc v 2 which corresponds to the standard Rayleigh- Jeans law. 

In addition to this, it is straightforward to generalize Stefan-Boltzmann law and show 
that it remains the same, i.e. it is still proportional to T 4 , but with a g-dependent constant 
{(T q ). To see this, let us write the total emitted power per unit surface, 

POO 

P q = / dvD q {y) . (20) 

J 

If we use here eq.(18) with the dimensionless variable x, 

p 9= ™$7*r — x -^- (2D 

h 3 c 3 Jo [1 - (1 - q ) x ]*=i -1 
can be obtained. Since the integral term is independent of T, it can be written down as 

P q = a q T A , (22) 



which is the generalized Stefan-Boltzmann law, with a g-dependent prefactor. 

It is also possible to obtain the generalized Wien shift law. By making eq.(18) maximum, 
v m (the frequency value which makes D q {y) maximum) yields a nonlinear equation : 



~kT 



(3g-4) + 3 



1 



hv m 



+ ^3(l-g)-3 = . 



(23) 



Although it is very difficult to find an analytical solution for this equation, the graphical 
solution is adequate for our purpose. From the graphical solution, for q = 0.95 and q = 1.05 
we have found hu m /kT = 2.444 and 3.347, respectively. For the same values of Tsallis 
g-index, eq.(19) in ref.[24] gives hu m /kT = 2.563 and 3.08. Once again it is possible to 
see that our results provide bounds to the exact ones. Moreover, the result v m {q > 1) > 
v m (q — 1) > v m {q < 1), which appears in ref.[24], is also valid in the present investigation. 
The important diversity between this result and that of the g^-deformed quantum groups 
investigations namely fm^gc) — V mi)-/Q.G) an d v m (q G < 1) < u m (q G = 1) is readily observed. 



IV. CONCLUSIONS 

Although the generalized distribution functions are first introduced in 1993, there has 
been no attempt to apply them to the physical systems until a very recent effort performed 
by Pennini et al. [39] where two single particles are considered. In this letter, our goal 
is to apply one of the generalized distribution functions (generalized Planck distribution) 
to another physical system (Planck law for the blackbody radiation). On the other hand, 
the approximate procedure used here to find a bound to the generalized Planck law is 
simpler than that of Tsallis et al. [24] (since the procedure is completely the same with 
that followed in standard textbooks of Statistical Physics for the derivation of standard 
Planck law) and also the bound seems to be more general, since not necessarily g = 1 (it 
will have a meaning, of course, if a physical system requiring this condition exists). Lastly, 
it is worthwhile to imply here that all the generalized laws derived here (i.e. generalized 
Planck law, generalized Rayleigh- Jeans law, Stefan-Boltzmann law and Wien law) transform 
to corresponding standard well-known laws in the q — > 1 limit. 
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FIGURE CAPTIONS 



Blackbody photon energy density per unit volume versus hu/kT in the frame of this work 
and ref. [24] . 
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Fig.l 




